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ABSTRACT 

The addition of Wide Field Camera 3 (WFC3) on the Hubble Space Telescope {HST) 
has led to a dramatic increase in our ability to study the z > 6 Universe. The im- 
provement in the near-infrared (NIR) sensitivity of WFC3 over previous instruments 
has enabled us to reach apparent magnitudes approaching 29 (AB). This allows us to 
probe the rest-frame ultraviolet (UV) continuum, redshifted into the NIR at z > 6. 
Taking advantage of the large optical depths of the intergalactic medium at this red- 
shift, resulting in the Lyman-a break, we use a combination of WFC3 imaging and 
pre-existing Advanced Camera for Surveys (ACS) imaging to search for z ?» 7 galaxies 
over 4 fields in and around Great Observatories Origins Survey (GOODS) South. Our 
analysis reveals 29 new z w 7 star forming galaxy candidates in addition to 15 pre- 
existing candidates already discovered in these fields. The improved statistics from our 
doubling of the robust sample of z-drop candidates confirms the previously observed 
evolution of the bright end of the luminosity function. 

Key words: galaxies: evolution - galaxies: formation - galaxies: starburst - galaxies: 
high-redshift - ultraviolet: galaxies 



1 INTRODUCTION 

There has been great progress in recent years in discover- 
ing star- forming galaxies at higli redshift [z > 5), either 
through tire Lyman-a emission line or their rest-frame UV 
continuum colours, which are strongly affected by absorp- 
tion by intervening hydrogen (the Gunn-Peterson effect at 
z > 6.3 provides almost complete absorption shortward of 
Lyman-Q, e.g. Gunn & Peterson 1965, Fan et al. 2001). This 
Lyman-break technique was initially used to identify z ~ 3 
galaxies from broad-band images (Steidel et al. 1996), and 
improvements in sensitivity of near-infrared imaging have 
enabled this to be pushed to higher redshifts (e.g. at z ~ 6 
using the Advanced Camera for Surveys on HST, Stanway, 
Bunker & McMahon 2003). 

Recently, the new Wide Field Camera 3 (WFC3) cam- 
era has been installed on HST, which has a near-infrared 
channel with significantly larger field and better sensitiv- 
ity than the previous-generation Near Infrared Camera and 
Multi-Object Spectrometer (NICMOS) instrument. A num- 
ber of multi-waveband imaging campaigns are underway, 
typically using broad-band filters at 1.0, 1.25 and 1.6/im 
and targetting fields with existing deep optical data. This 
is an ideal strategy to look for optical "drop-outs", seen 
only at the longer WFC3 wavelengths, which are candidate 
z > 7 galaxies. A number of papers (e.g. Oesch et al. 2010, 
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Bouwens et al. 2010, Bunker et al. 2010, McLure et al. 2010, 
Van et al. 2010 and Finkelstein et al. 2010, Wilkins et al. 
2010) have appeared in the past year presenting first results 
from such imaging. 

In this paper we analyse all the data taken so far in two 
large programs targetting fields in the vicinity of the Great 
Observatories Orgins Deep Survey-South (GOODS-South), 
which in one case includes the Hubble Ultra Deep Field, and 
present a new sample of z-band drop-outs which more than 
doubles the existing number of robust z > 7 candidates. By 
building such large samples (> 40 objects spanning 3 magni- 
tudes) we can ultimately study the rest-UV luminosity func- 
tion, and hence estimate the integrated star formation rate 
within 800 Myr of the Big Bang. Understanding the global 
star formation history at these redshifts is crucial in answer- 
ing the question whether the UV photons produced by the 
short-lived massive OB stars were sufficient to reionize the 
Universe. 

This paper is organised as follows: in Section 2 we dis- 
cuss the observations, our data reduction process and our 
photometry measurement and catalogue production tech- 
nique. In Section 3 our candidate selection procedure is de- 
scribed and we present a list of z-drop candidates together 
with a comparison with other studies. In Section 4 we in- 
vestigate the properties of these objects, including their red- 
shift distribution, surface density, and luminosity function. 
Further, in this section we investigate the implications of 
the luminosity function for the star formation rate density 
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and production of ionising photons. In Section 5 we present 
our conclusions. Throughout, we adopt the standard concor- 
dance cosmology of Qm — 0.3, Qa = 0.7 and use Ho = 70 
kms"'^ Mpc~^. All magnitudes are on the AB system (Oke 
& Gunn 1983). 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Observations 

In this paper we analyse ACS 6435™ , weoe™ , 1775™ and zs50ip- 
band imaging together with recently obtained near-infrared 
images from WFC3 on HST in Yo9Sm/Yio5w, Ji25w and 
J^i6o™-bands (see Figure 1, 1-1.6/im) taken in the vicinity 
of CDFS/GOODS-South. 

The WFC3 data come from two different HST 
programs. The Early Release Science (ERS) program 
GO/DD#11359 (RI. R. O'ConneU) covers ten overlapping 
pointings with two orbits in each filter; analysis of the first 
6 pointings was presented in Wilkins et al. (2010), and we 
here we analyse the complete dataset (referred to as the 
ERS field). Secondly, deep WFC3 imaging of the Hubble 
Ultra Deep Field (HUDF, Beckwith et al. 2006) and two 
nearby deep flanking fields come from the HST Treasury 
programme GO-11563 (P.I. G. Illingworth) . These flanking 
flelds were imaged by ACS in veoew, irrsm- and zgsoip-bands 
during 2005-6 in parallel with deep HST NICMOS NIC3 
observations of the original UDF as part of program GO- 
10632 (P.I. M. Stiavelli), and we adopt the original nam- 
ing convention from the UDF parallels of HUDF (NIC)P12 
and HUDF (NIC)P34 (as in Oesch et al. 2007, 2009), which 
we abbreviate to P12 and P34. We note that P12 is called 
HUDF05-01 in the WFC3 programme GO-11563, and 
P34 is HUDF05-02. The three WFC3 flelds observed faU en- 
tirely within the 11 arcmin? footprint of the ACS HUDF 
and flanking field images. This 192-orbit WFC3 program is 
ongoing, with the HUDF observations comprising half of the 
total orbits. In Bunker et al. (2010) we analysed the HUDF 
data obtained soon-after the commissioning of WFC3, and 
in this paper we study the new images of the two deep flank- 
ing flelds, and reanalyse the HUDF data using more recent 
on-orbit calibration of the WFC3 IR detector. 

The infrared channel of WFC3 was used, which is a 
Teledyne 1014 x 1014 pixel HgCdTe detector, with a field 
of view of 123" x 136" (a 10-pixel strip on the edge is 
not illuminated by sky and used for pedestal estimation). 
Both programs use the same J and //-band filters (F125W 
and F160W), although the ERS images use a F-band filter 
(F098M) which covers only the blue side of the wider F105W 
filter used in the HUDF and fianking field images (see Figure 
1). The data were taken in "MULTIACCUM" mode using 
SPARSAMPLEIOO, which non-destructively reads the array 
every 100 seconds. These repeated non-destructive reads of 
the infrared array allow gradient-fitting to obtain the count 
rate ("sampling up the ramp") and the flagging and rejec- 
tion of cosmic ray strikes. 

For the HUDF and parallel fields, each MULTIACCUM 
exposure comprised 16 reads for a total duration of 1403 sec 
per exposure, with two exposures per orbit. The WFC3 
HUDF data, taken over 26 August - 05 September 2009 
U.T., comprises 39.2 ksec (20 exposures) in Yiosm (not in- 
cluding 8 exposures affected by severe persistence. Bunker 



et al. 2010), 44.8 ksec (32 exposures) in Ji25m and 78.5 ksec 
(54 exposures) in Hieow- The new P34 data (spanning 23 
January - 17 February 2010 U.T.) goes almost as deep, 
with 20, 28 & 34 exposures in Yiosiu, J1251U & Hieow respec- 
tively (28.0 ksec, 39.2 ksec & 47.7 ksec); observations of this 
fields are nearly complete, unlike P12 where the bulk of the 
data will come next year. In this paper we analyse the exist- 
ing 12 exposures (16.8 ksec) in Yiosui-band and 24 exposures 
(33.6 ksec) in Ji25iu-band for P12 (taken over 02 - 15 Novem- 
ber 2009 U.T.), although only 4 exposures (5.6 ksec) have 
been taken so far in the //leom-band. For the ERS data, each 
MULTIACCUM exposure comprised 9 or 10 non-destructive 
read-outs, totalling 803-903 s, and there were three such ex- 
posures per orbit (and each filter was imaged for 2 orbits 
per pointing, with the 6 exposures totalling 2.6 ksec). 

Advanced Camera for Surveys (ACS) optical/NIR 
(0.4 — l.Ofim) imaging for the 4 fields comes from a variety 
of sources. For the ERS field we utilise 6435^,, veoe-w, 1775™, 
and zgsoip-band imaging from GOODSv2.0 (Giavalisco et al. 
2004) . ACS 6435™ , V606w , irTBn, , and zgsoip-band Imaging of 
the HUDF is taken from Beckwith et al. 2006. For the P12 
field we use the publicly-available vqo6w , irrsw , zgsoip reduc- 
tions provided by the UDF05 team (Oesch et al. 2007) and 
for the P34 field we obtain veosm , irrsm , ^gsoip-band imaging 
from the HST archive and reduce it as described below. 

2.2 Data Reduction 

The WFC3 imaging was reduced by first utilising the 
IRAF.STSDAS pipeline calwf c3 to calculate the count rate 
and reject cosmic rays through gradient fitting, as well as 
subtracting the zeroth read and fiat-flelding. We used MUL- 
TIDRIZZLE (Koekemoer et al. 2002) to combine exposures 
taken through the same filter in each pointing, taking ac- 
count of the geometric distortions and mapping on to an 
output pixel size of 0'.'06 from an original 0'.'13pix~^. This 
was the same scale as we used in our analysis of the Hubble 
Ultra Deep Field WFC3 images (Bunker et al. 2010) and 
corresponds to a 2 x 2 block-averaging of the GOODSv2.0 
ACS drizzled images and HUDF ACS reduced images. 

The new data were reduced using the latest version 
of calwf c3 (29 October 2009 release), and we also re- 
reduced the HUDF images originally presented in Bunker 
et al. (2010) - the earlier paper had used a reduction us- 
ing XDIMSUM and fits to the geometric distortion due to 
MULTIDRIZZLE not then being available for the newly- 
commissioned WFC3. For the HUDF and fianking fields, 
we used a DRIZZLE pixel fraction of 0.6 to recover some 
of the under-sampling. In each of these three fields we sur- 
vey 4.18arcmin^ in all exposures, with another 0.67arcmin^ 
surveyed at half the maximum depth. For the ERS field, we 
reduced the data for all 10 ERS pointings using the same 
technique. However, we used a pixel fraction of 1.0 in mul- 
tidrizzle, as we only had 6 exposures. We then mosaicked 
together the 10 ERS pointings in each fllter, using inverse- 
variance weighting for the overlap regions, producing a held 
of fairly uniform depth covering 37arcmin^, with a further 
8arcmin^ going less deep. 

ACS V606-W , i775w , zgsoip-band imaging of the P34 held 
taken from the HST archive was reduced using MUL- 
TIDRIZZLE to combine blocks of data taken close in time 
with similar telescope roll angles, using an output 0'.'03 pixel 



New Star Forming Galaxies at z ^ 7 3 



ACS/WFC3 2(T Detection Limit 



Field ID 


Centre (J2000) 


Area 


f'435m 


V606w 


«775tD 




^098m/105m" 


Jl25w (7o-) 




HUDF 


03:32:38.4 -27:47:00 


4.2 


30.3 


30.7 


30.6 


30.0 


29.65 


29.72 (28.36) 


29.67 


P34 


3:33:05.3 -27:51:23 


4.2 




29.9 


29.5 


29.7 


29.45 


29.58 (28.22) 


29.41 


P12 


3:33:01.9 -27:41:10 


4.2 




29.9 


29.6 


29.6 


29.16 


29.50 (28.14) 


28.23 


ERS 


3:32:23.6 -27:42:50 


37.0 


29.1 


29.1 


28.5 


28.4 


28.02 


28.39 (27.03) 


28.10 



Table 1. Summary of observations. "VogSm for the ERS field and Yios™ for the HUDF/P12/P34 fields measured in 0'.'6 diameter 
apertures. 



scale. These subsets of drizzled images were then registered 
and combined with IRAF . imcombine. Our combined images 
were 4.8ksec in u-band, 10.6 ksec in i-band and 26.8 ksec in 
z-band. All the ACS images were then block-averaged 2x2 
and registered with our drizzled WFC3 frames. 

2.2.1 Photometry and Zeropoints 

In our final combined Vogsm/iosiD-band images, we measure 
a FWHM of ~ 0'.'15 for point sources in the field. As most 
high-redshift galaxies are likely to be barely resolved (e.g., 
Bunker et al. 2004, Fergusson et al. 2004) we perform pho- 
tometry using fixed apertures of 0'.'6 diameter, and intro- 
duce an aperture correction to account for the flux falling 
outside of the aperture. This correction was determined to 
be 0.2 mag in the yo98m/i25m-band and 0.25 mag in J12511,- 
band and ffieoiu-bands from photometry with larger aper- 
tures on bright but unsaturated point sources. We note that 
the 7ifi60m-band images display significant Airy diffraction 
rings around point sources. For the ACS images, the better 
resolution and finer pixel sampling require a smaller aper- 
ture correction of « 0.1 mag. 

For the ACS images we take the standard ACS zero- 
points. For WFC3, we use the recent zeropoints reported 
on 'http : / /www . stsci . edu/hst7wf c3/phot_zp_lbn during 
February 2010, where the Yogsm-band has an AB magni- 
tude zeropoint of 25.68 (such that a source of this bright- 
ness would have a count rate of 1 electron per second), and 
zeropoints of 26.27, 26.25 & 25.96 for Yios™, Ji25m & Hiaow. 
We note that the information in the image headers of the 
earlier images (released in September 2009) is slightly dif- 
ferent by 0.1-0.15 mag, with zeropoints of Viosiu = 26.16, 
Ji2£,w = 26.10 & Hieom = 25.81 (as used in Bunker et al. 
2010). 

We also apply a correction for the small amount of fore- 
ground Galactic extinction toward the GOODS-South using 
the COBE/DIRBE & IRAS/ISSA dust maps of Schlegel, 
Finkbeiner & Davis (1998). The optical reddening is E{B — 
V) = 0.009, equivalent to extinctions of Agsoip = 0.012, 
^105™ = 0.010, = 0.008 & Ai6o» = 0.005. 

Information, including the position, area, and average 
2a depths reached in each filter, for each field are given in 
Table 1. 

2.3 Construction of Catalogs 

Candidate selection for all objects in the field was performed 
using version 2.5.0 of the SExtractor photometry package 
(Berlin & Arnouts 1996). For the zssoip-drops, as we are 
searching specifically for objects which are securely detected 
in the WFC3 Yo98m/io5uj-band, with minimal flux in the 
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50 Uyr Cont. SF. Z=Z^. z=8.0 




3.6 3.8 4 4.2 

log,Q(X[A]) (observed frame) 



Figure 1. Top panel - Model (from the Starburst99, Leitherer et 
al. 1999) speetral energy distribution (SED) of a redshifted 2 = 7 
star forming galaxy. Middle panel - Potential contaminants: Ob- 
served SED of a low-mass dwarf star (class: r4.5, Knapp et al. 
2004) together with the model (Starburst99) SED of a 3.5Gyr 
Single-aged Stellar Population at z = 1.8. The bottom two pan- 
els show the transmission functions of the combination of filters 
available to each field. 

ACS images (see following section), fixed circular apertures 
0'.'6 in diameter were 'trained' in the Voggm/iosm-image and 
the identified apertures used to measure the flux at the same 
spatial location in the remaining ACS and WFC3 bands by 
running SExtractor in dual-image mode. For the initial ob- 
ject identification, we adopted a limit of at least 5 contiguous 
pixels above a threshold of 2a per pixel (on the data drizzled 
to a scale of 0'.'06 pixel" ). 



3 CANDIDATE SELECTION 

At very high redshifts the integrated optical depth of neu- 
tral hydrogen in the intergalactic medium (IGM) becomes 
large, resulting in a Gunn-Peterson trough (Gunn & Peter- 
son 1965), a large flux decrement shortward of the rest-frame 
wavelength of the Lyman-a transition (1216A). At z ~ 7 the 
location of the Lyman-a break is redshifted to ^ l^m - the 
ACS Zssoip and WFC3 Yio5u./0987n are suitably located such 
that a, z = 6.5 — 8.0 star forming galaxy will experience a 
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significant flux decrement between the two filters (see Fig- 
ure SpecificaUy a jzgsoip — Yiosw > 1.0 colour is indicative 
of such a spectral break. 

In order to select only robust 2 ~ 7 star forming galaxy 
candidates we first introduce a J12511, magnitude cut. The 
choice of the J12511, filter for this cut is motivated by the fact 
that for z ^ 7 galaxies this filter will be uncontaminated 
by either the break or Lyman-a emission and will thus be 
an indication of the rest-frame ultraviolet luminosity of the 
candidates. For each field we define the J12511, limit as the 
magnitude where the mean signal-to-noise ratio of an object 
at the limit is 7a. 



3.1 Contamination 

When searching for distant galaxies using only broadband 
photometry, contamination becomes a serious issue. Sources 
of contamination can be placed into three broad categories: 
objects whose intrinsic colours are similar to those of the 
target population; objects whose observed colours are simi- 
lar because of photometric noise; and the effect of transient 
phenomena. We discuss each of these below. 



3.1.1 Intrinsically Red Objects 

Two distinct types of contaminating objects can have ap- 



parent zssoip - 



105u)/098 



colours consistent with those 



of z ~ 7 galaxies, making them indistinguishable with- 
out the aid of additional constraints. These include lower- 
redshift galaxies where the zg^oip and yiosm/ogsm filters 
straddle the Balmer/4000A break feature (at approximately 
z 1.5 — 1.9) and low-mass dwarf stars, particularly those of 
the L and T classification whose low temperature and broad 
absorption features can mimic a spectral break. 

Examples of the spectral energy distributions (SEDs) 
of each of these types of object (a model 3.5 Gyr old single- 
aged stellar population at z = 1.8 and a T4.5 dwarf star) are 
shown in Figure [T] In both these cases the slope of the SED 
longward of the spectral break (i.e. longward of zgsoip) is 
somewhat redder than that predicted for a high- 2: star form- 
ing galaxy. The addition of a further filter at wavelengths 



redder than the Y^ 



105™/098m 



filter can then be used to dis- 



forming galaxies are shown in the (^gsoip — Yj 

w 1 098m Ji; 



criminate between high-2 galaxies and these lower-redshift 
interlopers. 

The ability for an additional redward filter to discrim- 
inate against these interlopers can be seen more clearly 
in Figures [2] and [S] In these figures the positions of both 
the interlopers and the tracks expected for high-2 star 

l()5tu/098m) ~ 

colour plane. With the exception of 
the lowest temperature T dwarfs where the Voggm filter is 
employed (the ERS field) these interlopers form a distinct 
locus separate from z k, 7 star forming galaxies with UV 
spectral slope index /? < 0.0 (where j\ — \^ is used as a 
model of the UV properties of star forming galaxies). 

Using this analysis it is possible to design a window in 

(285Wp - i^l05™/098m) " (^105^/098™ " ^125™) Colour - Colour 

space that selects only high-z star forming galaxies. For the 
HUDF/P12/P34 fields (i.e. where we have yios™ imaging) 
this zY J selection selection criteria is: 

{zaUMp - VioSuj) > 1.0 




(^105 '.'i25)a6 





Figure 2. liosiu — Jl25m and 2850;p — ^lOSii. colour - colour figures 
for the HUDF/P34/P12 fields showing our zY J colour selection 
window (grey shaded area), the location of our candidates, the 
predicted paths taken by high-2 galaxies (blue lines, /3 = —3.0 
(left) and /3 = 0.0 (right)) and the location of possible contami- 
nating sources. Contaminating sources include galactic stars (de- 
noted by coloured filled circles) and a passively evolving instanta- 
neous burst of star formation that occured at z = 10 (red dashed 
line). High-2 candidates arc denoted by black circles (where the 
size of the circle is an indication of the apparent Ji25w magni- 
tude). Small blue circles denote objects that meet our zYJ colour 
criteria but which are detected in the bluer (b)vi ACS bands. 
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■ Observed: J,2s,<26.0 
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■ SimuioLed: J,25,<i^b-D 


■ '7 ■ 

■ 'i i'i .' i - 


■ 1 ■ 



-0.5 0.5 -0.5 0.5 

(^105~^125)aB 

Figure 4. Comparison of tlie observed (left panel) (zgsoip — 
YiOSw) - (^lOSm — •^I25m) colours of objects in the HUDF with 
that simulated by adding noise to observed distribution of bright 
galaxies (right panel). In both cases objects meeting our zYJ 
criteria are indicated by the enlarged points. 



Figure 3. Yo9Sm — Ji25w and zssoip— ^09gm colour - colour figures 
for the ERS field showing our zYJ colour selection window (grey 
shaded area), the location of our candidates, the predicted paths 
taken by high-z galaxies (blue lines, /3 = —3.0 (left) and /? = 0.0 
(right)) and the location of possible contaminating sources. Con- 
taminating sources include galactic stars (denoted by coloured 
filled circles) and a passively evolving instantaneous burst of star 
formation that occured at 2: = 10 (red dashed line). High-2 can- 
didates are denoted by black circles (where the size of the circle is 
an indication of the apparent Ji25uj magnitude). Small blue cir- 
cles denote objects that meet our zYJ colour criteria but which 
are detected in the bluer bvi ACS bands. 



{zsrMp - YuKw) > 2.4 X (Fiosm - J125™) + 0.4 

(VlOSuj — Jl25iu) < 1-0 



The use of an alternative Y filter (Vogsm) in the ERS 
field necessitates the use of a slightly different criteria: 

{zS50lp — ^098m) > 1-2 

{zS50lp — ^098m) > 0.9 X (io987n — Jl25w) + 0.7 
(i^098m — Jl25w) < 2.0 

It is important to note that such criteria biases against 
a hypothetical population of z ~ 7 galaxies which have ex- 
tremely red UV spectral slopes (resulting in Ji25m — -ffieou; > 
1.0 colours), which could be attributed to the presence of 
large amounts of dust attenuation. The potential bias of 
our selection criteria and a more general analysis of the UV 
properties of the candidates presented in this work is dis- 
cussed in more detail in Wilkins et al. (2010b) where we 
conclude that the distribution of UV spectral slope indices 
is consistent with /3 « —2.2 (i.e. rather blue). 

A second consideration to note is that while our selec- 
tion criteria excludes all known intrinsic contaminants, it 
remains possible that there exists a hitherto unknown pop- 
ulation of objects whose intrinsic photometry may be con- 
sistent with our selection criteria. 



3.1.2 Photometric Scatter 

While the selection criteria described in the preceding sec- 
tion is able to guard against contaminants whose intrinsic 
photometry matches that which is observed, photometric 
noise may scatter objects into our selection window, provid- 
ing a significant source of contamination. 

To asses the impact of photometric scatter when work- 
ing at low signal-to-noise, we performed a simple simulation 
to generate predictions at fainter magnitudes based on the 
colour distribution of more luminous objects. Using the ob- 
served bvizYJH photometry distribution of bright sources 
(24.0 < Ji25» < 26.0 in the case of the HUDF) synthetic 
objects were inserted into the original images with fainter 
magnitudes. The number of inserted synthetic sources in a 
given magnitude interval was chosen to recreate the observed 
number of sources in the same interval. 

Selecting objects in the HUDF with Ji25u, < 28.5, sug- 
gests, from this simulation, that ^ 14 objects would scatter 
into our zY J selection window, surpassing the number of po- 
tential high-z candidates. This is shown in Figure |4] where, 
in the right panel, the simulated colours of Ji25uj < 28.5 
objects are shown and in the left the observed distribution 
(again Ji25m < 28.5) 

To guard against this serious level of contamination we 
take advantage of the deep optical imaging that is avail- 
able in the ACS 6435™, ueoe™ and 177511, bands. Star forming 
galaxies at z ~ 7 axe not expected to contribute any flux in 
these bands as the lie blueward of the rest-frame location of 
the Lyman-a 1216A, thus we impose an additional bvi non- 
detection criteria for the selection of our candidates. That 
is, we classify any object with a > 2(j detection in any of 
the 6435^ , vqo6w and 177511, filters as a lower-redshift contam- 
inant. Imposing this selection criteria into our simulation we 
then find significantly fewer potenitial contaminants (~ 0.6 
in the HUDF). 

Making the J125U, magnitude limit fainter by 0.5 mag 
results in the contamination rate jumping significantly. This 
analysis highlights the necessity of having additional deep 
optical imaging to remove any such contamination. In total, 
over all fields, we expect a contamination rate of ~ 2.4 ob- 
jects due to photometric scatter (or a ~ 5% contamination 
rate) (assuming a Ji25uj limit corresponding to the average 
6(T depth). 
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Figure 5. 2'.'. 4 X 2'.'. 4 composite (B=v, G=iz, R=YJ) and izY 
thumbnail images of a significant apparent motion object in the 
P12 field. 



3.1.3 Transients 

The large time delay between the acquisition of the ACS 
zssoip (taken between 2004 and 2006) and the WFC3 YJH 
imaging leaves open the possibility of contamination by tran- 
sient phenomena such as supernovae (SNe) and objects with 
significant apparent motions. Indeed, it appears both these 
phenomena are contaminants within our fields. 

Supernovae occurring within the WFC3 observation 
window could produce significant flux in the YJH bands. 
This would naturally produce a large flux increase between 
the ACS bviz-haiids (in which the transient would not con- 
tribute) and the NIR filters, mimicking the effect of a (poten- 
tially) very strong spectral break. The analysis of Bouwens 
et al. (2008) suggests a supernovae contamination rate of 
0.012 per arcmin^ suggesting that we would expect to see 
~ 0.75 SN over aU the fields. 

Contamination from supernovae can, however, at the 
bright end be guarded against by comparing the current 
NIR imaging with pre-exist NIR imaging. The brightest pre- 
liminary candidate in the HUDF (which is not included in 
our final list of candidates), with Ji25w = 25.6 (fiagged as 
zDO in Bunker et al. 2010 and noted in Oesch et al. 2010, 
McLure et al. 2010, Finkelstein et al. 2010 and Yan et al. 
2010) should be detected in previous NICMOS imaging of 
the field (which has a 2a depth of Jnom ~ 26.7) had it been 
a persistent object (i.e. not a transient). Its absence in this 
additional imaging is used as evidence that it is a contam- 
inant and it is thus excluded in our analysis (as in Bunker 
et al. 2010). 

Phenomena with significant apparent motions are pos- 
sible contaminants as the position of the object when the 
z-band imaging was acquired will be offset from that when 
the WFC3 NIR imaging was obtained. Photometry based on 
the location of the source in the Y-band image could then ar- 
tificially suggest a break between the WFC3 and ylCS* bands. 
Through visual inspection of the zssoip and VogsTn/ios™ i™" 
ages for each candidate such objects can be guarded against. 
One such object was found in the P12 field (thumbnails of 
the RGB composite and izY images is shown in Figure [5]) 
and was removed from subsequent analysis. 



3.2 



7 Candidate Galaxies 



Using the selection criteria described in the preceding sub- 
section we compile a list of candidate z ~ 7 star forming 
galaxies in the HUDF, P34, P12 and ERS fields. The po- 
sitions and photometry of these objects are listed in Table 
[2] Thumbnails of the bvizYJH images of these candidates 
(where available) are presented in Figures [§]- [HI In total we 
find 44 (HUDFTl, P34:15, P12:7, ERS:11) candidates cov- 
ering a range of apparent Ji25id magnitudes of 25.9 — 28.36. 
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Figure 6. 2'.'4 X 2''4 bvizYJH thumbnail images of objects meet- 
ing our selection criteria in the HUDF field. 



A cross-matching of these candidates with those presented 
in other studies is undertaken in Section 4. 



3.3 Comparison with Other Studies 

Both the WFC3 HUDF data (Bunker et al. 2010, Finkelstein 
et al. 2010, McLure et al. 2010, Oesch et al. 2010 and Yan 
et al. 2010) summarised in Table (H)) and some of the ERS 
pointings (Wilkins et al. 2010) have been previously anal- 
ysed. In each case different reduction techniques, methods 
for measuring the photometry (including the choice of pho- 
tometric zeropoints) and 2 f» 7 selection criteria were used. 
A brief summary of these studies is given in Table [3] where 
information about their selection methodology and criteria 
is documented. 

In the case of ERS dataset (partially analysed by 
Wilkins et al. 2010) we match 4 of the 6 candidates (1, 2, 
4 and 6) with the remaining two failing to make our Ji25m 
selection criteria (Wilkins et al. 2010 employed a Yogsm cri- 
teria). 

Cross-matching our candidates with those presented in 
the studies listed in Table [3] indicates that all our HUDF 
candidates exist in at least two of the other studies. A cross 
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Table 2. z -band drop out candidate z K, 7 galaxies meeting our selection criteria in the HUDF, P34, P12 and ERS fields. Objects 
are divided by field and then ordered by apparent Ji25m magnitude. Photometry based on 0'.'.6 diameter apertures with an aperture 
correction. Quoted limits correspond to 2a. " We introduce a small offset to the astrometry in the HST image headers of P34 of 
ARA = 0!03 to match the astrometry from GSC-2 and 2MASS. ^ This object has a J125 magnitude 0.1 fainter than the adopted 7cr 
limit and was erroneously included in the initial draft of this work. It is omitted from all subsequent analysis. 



Study 



Method 



magnitude limit 



This Study 
BIO" Bunker et al. 2010 
OlO" Oesch et al. 2010 
MIO" McLure et al. 2010 
FIO Finkelstein et al. 2010 
YIO Yan et al. 2010 



z — Y+Y — J colour criteria + hvi < 2(t 
z — Y+Y — J colour criteria + bvi < 2(t 
z — Y+Y — J colour criteria + bvi < 2cr 
bvizYJH photometric redshifts 
bvizYJH photometric redshifts 
z — Y colour criteria bvi < 2cr 



Ji25m < 28.5 
Yiost^, < 28.5 

S/N{Yi05^ A Jl25ru) > 5o- 
S/N{Yl05w V Jl25m V HlliOm) > 5a 

S/N{Ji25w A Hieow) > 3.5o- 
S/N{Y-io5w) > 3.0 



Table 3. Summary of other studies that analysed the HUDF field. " Employed different zeropoints. 
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40566437 
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Table 4. A cross-matched Hst of our candidates and those reported by the studies listed in Table [3] and any additional studies. Numbers 
in parentheses in the McLure et al. (2010) (MIO) and Finkelstein et al. (2010) (FIO) columns denote the photometric redshifts recovered 
by these studies. 



P34.Z.2397 I 




P34.Z.6 



P34.Z.3900 ] 
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« 



Figure 7. 2'.'4x2'.'4 vizY JH tliumbuail images of objects meeting 
our selection criteria in the P34 field. 




Figure 8. 2'.'4 X 2'.'4 vizYJH thumbnail images of objects meeting 
our selection criteria in the P12 field. 



matched list of our candidates and the corresponding iden- 
tifications from the other studies are shown in Table |3] 

Due to differences in the selction technique (and cri- 
teria), there are several candidate objects included in the 
studies listed in Table [3] which fail to enter our selection. 
These are listed and cross- matched in Table [S] The McLure 
et al. (2010) and Finkelstein et al. (2010) analysis used a 
photometric redshift selection technique; as such their cat- 
alogues may include objects at too high, or too low redshift 
to meet our selection criteria, depending on their choice of 
redshift range. In Tables |4] and [5] only objects with photmet- 
ric redshifts in the range 6.3 < z < 7.8 are included. This 
range is motivated by a series of simulations discussed in 
more depth in Section [4.21 

A significant fraction of the objects included in Table [5] 
are excluded from our selection on the basis of our conser- 
vative Ji25ijj limit. For example two objects (zDlO and zD7) 
in Bunker et al. (2010) (who adopt a Viosii. magnitude limit 
criteria) are excluded on this basis (5/8 remaining objects 
from Oesch et al. 2010 are excluded for this reason, as are 
4/11 of the remaining Yan et al. 2010 objects). 
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ERS.Z.4022S 
ERS.Z.20851 




Figure 9. 2'.'4x2''4 bvizYJH thumbnail images of objects meet- 
ing our selection criteria in the ERS field. 



The remainder of the excluded objects fail to meet our 
izs50ip - Yws-w) ■■ {Yl05^u - Ji25m) colour Selection criteria 
(i.e. none are excluded on the basis of an optical detection). 
Bunker et al. (2010) initially employed a zssoip — Viostu cut 
before considering the Viosm — J1251U colour to exclude very- 
red potential low-redshift and stellar interlopers. As a result 
several objects (zDll, zD9, zD4) which make our Ji25-w se- 
lection criteria fail to make our colour selection criteria. On 
inspection of the photometry of these objects they are re- 
vealed to all lie very close to our selection window suggesting 
they are potential high-redshift galaxies and not contami- 
nants. 

The exlcusion of these candidates can be seen more 
clearly in Figure 1101 In this figure the points denote the lo- 
cation (in the zssoip - Yiosw - Yiosn, - Ji25w colour plane) of 
candidates from the Oesch et al. (2010), Bunker et al. (2010) 
and Yan et al. (2010) studies and the corresponding arrows 
point to the location based on our photometry. Those ob- 
jects meeting our selection criteria based on our photometry 
are denoted with a halo. The figure shows that candidates 
failing to make our selection criteria generally lie close to the 
selection window. Figure [10] also highlights that there exist 
sometimes significant differences between the photometry of 
the various studies. To some extent this is due to differences 
in the assumed photometric zeropoints (i.e. no attempt was 
made to harmonise these) with the remainder of the devia- 



>- 
I 



T-r-r-T-p-r-r-r-r-p-r-r 

.Bunker et al. 2010 



I I I I I I I I I I I I I I 

_ian et ol. 2010 




.Oesch et al. 2010 




Figure 10. Comparison of the zgsoip — ^105™ and Yiosii. — Ji25w 
colour for candidates reported by other studies from both the 
original photometry (i.e. from the other studies) and our pho- 
tometry. In each case the original photometry is denoted by the 
position of the circle. The associated arrows point to the position 
in colour-colour space where the object would lie according to 
our photometry. Objects included as candidates in this study are 
circled. 



tions due to differences in the reduction process and choice 
of photometry methodology. 

Comparing our candidate list with those of McLure et 
al. (2010) and Finkelstein et al. (2010) is slightly more com- 
plicated as they employ photometric redshift techniques and 
present candidates encompassing a larger redshift range. Of 
their candidates lying within our predicted redshift selection 
window (6.3 < z < 7.8 , see Section [4. 2p several have Ji25m 
magnitudes fainter than our imposed limit. The remaining 
objects have {zssoip - Viosii.) : (Vios™ - Ji25-w) colours which 
are typically only marginally inconsistent with our selection 
criteria. 

We have imposed a colour selection criteria and S/N 
limit where contamination from photometric scatter is min- 
imal. Other groups have presented objects outside our colour 
selection, or at lower S/N which may plausibly be high- 
redshift, but have a greater possibility of being contami- 
nants. The omission of these candidates from our selection, 
however, does not impact any of the conclusions drawn from 
the observed luminosity function (see Section I4.3|l . As de- 
scribed in the following sections, the luminosity function is 
determined by considering the effective volume of the survey, 
which is sensitive to the choice of selection criteria. Thus, 
our conservative selection criteria should recover the true 
luminosity function of 2; « 7 star forming galaxies. 



4 CANDIDATE PROPERTIES 
4.1 Surface Densities 

One useful diagnostic of the z-drop population probed in 
each field are the observed surface densities as a function of 
limiting Ji25uj-magnitude (as shown in figure [11] for each 
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Table 5. List of candidates from the studies outlined in Table |3] that fail to make it into our selection together with the reason why 
(final column). A - fails to meet our Ji25u,-band magnitude cut. B - fails to meet our colour selection, f would have met our selection 
criteria if photometry from original study was adopted. 

potential of high-redshift z-drop galaxies as contaminants 
in, for example, the search for very low-mass stars (e.g. De- 
lorme et al. 2010). 

The observed surface densities are consistent between 
fields (i.e. within the uncertainties) with the exception that 
the P34 surface dnesities exceed those of the ERS. How- 
ever, this is unsurprising as the surface density of z-drops 
in the ERS is not directly comparable to that the other 
fields due to the use of the Fo98m filter (over the Fios™ used 
in the HUDF/P34/P12 fields). The agreement between the 
3 deeper fields (HUDF/P34/P12), which are seperated by 
~ 2 Mpc suggests we are not probing a compact cluster 
which would affect our estimates of the rest-frame UV lu- 
minosity function and thus star formation rate density and 
ionising photon luminosity density. 

4.2 Completeness and Redshift Distributions 

A useful description of our z-drop population is the theo- 
retical redshift distribution we would expect to probe, that 
is, it is useful to ascertain the probability of recovering a 
high-2 galaxy as a function of redshift and rest-frame UV 
luminosity - p(Mi6oo, z). 

To quantify this probability a series of simulations is 
conducted for each field. These simulations are performed by 
inserting sources with properties similar to those expected 
for high-z galaxies (i.e. compact sources with half-light radii 
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Figure 11. The observed cumulative surface density (i.e. the 
number of objects brighter than the limiting Ji25m of z « 7 can- 
didate galaxies for each of our fields. Note, because of the differ- 
ent Y filters employed by the HUDF/P34/P12 and ERS fields the 
surface densities are not strictly comparable due to the slightly 
different redshift selection. 



field). This gives a direct indication of any field to field 
variance due to cosmic variance, that is whether one field 
is probing an area of higher-density than the average - a 
potential problem for narrow pencil-beam surveys. Further, 
knowledge of surface density is useful when considering the 
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Figure 12. The probability of recovering a galaxy in our simu- 
lation for each field as a function of redshift for several different 
absolute rest-frame Mieoo magnitudes. The mean redshift is in 
each case denoted by a dot. 



~ 0.1") into the original images and running our selection 
process to determine the probability that a source would be 
recovered (and included in our selection). 

A further important parameter in the simulation is 
the choice of the UV continuum properties. This is pa- 
rameterised throught the UV spectral slope index /3 (where 
/a cx: over the rest-frame UV continuum). Due to our 
selection criteria, objects with different values of /? enter the 
selection window at slightly different redshifts which would 
result in a different p{Mi(joo,z). However for objects with 
—3.0 < /3 < 0.0 the redshift at which objects enter the 
selection window is fairly uniform and the resulting values 
of p{MiQoo,z) are very similar. Objects with /3 > 0.0 (sug- 
gestive of large dust attenuations) are however excluded as 
part of our strategy to eliminate low redshift contaminants. 
However previous studies (e.g. Bouwens et al. 2009) have 
found that the UV slopes of high-redshift objects are typ- 
ically blue with /3 ~ —2.2, corresponding to a star form- 
ing population which is essentially dust-free (under common 
assumptions regarding the previous star formation history, 
metalhcity and initial mass function). For our candidates we 
find a mean {Ji2£,w — -ffieom) colour of 0.0 corresponding to 
13 = —2.0 and find few objects with (J125™ — Hieom) > 0.25 
(i.e. /3 > —1.0). In our simulations we assume that the UV 
spectral slopes of objects are distributed as a gaussian with 

= -2.2 and a = 0.5. 

The results of these simulations are shown in Figure [121 
In the HUDF and flanking fields galaxies with large intrin- 
sic luminosities (Mieoo < —20.0) p{Mieoo,z) remains fairly 
constant over 6.3 < z < 7.8. At lower luminosities the prob- 
ability of recovering a galaxy quickly decreases at the high- 
z end of this range due to increasing luminosity distance 
and Lyman-alpha forest blanketting of the Fios^-filter with 
redshift. The ERS field, by virtue of the narrower VcQSm , 
has a slightly more restricted redshift probability function 
than the other fields. Irrespective of the field, though with 
a slight bias to a higher-redshift for intrinsically more lu- 



minous galaxies, the expected mean redshift of zssoip-drop 
high- 2 star forming galaxies is z ~ 7.0. 



4.3 Luminosity Function 

The probability of observing a galaxy as a function of ab- 
solute rest-frame UV magnitude and redshift, ^(Migoo,^), 
described in the preceding section, can be used to generate 
a prediction of the effective volume surveyed as a function 
of absolute Mieoo magnitude which can then be used to esti- 
mate the piecewise luminosity function (LF). The relation- 
ship between the p{M, z) and effective volume Vcff is given 
by 



p{M,z)—dz, 
dz 



(1) 



as described in (e.g. Steidel et al. 1999, Stanway et al. 2003). 

To estimate the LF we assume that all galaxies lie at the 
centre of the computed redshift distribution and determine 
the absolute magnitude. At z = 7.0 the effective wavelength 
of the Ji25ii. filter corresponds to a rest-frame value of 1563A. 
For ease of comparison with other studies we fc-correct the 
observed J1251U magnitudes to I6OOA by using the observed 
(Ji25iu — HiQOw) colour (though it is worth noting that this 
is a small correction given that (Ji25m — Hi^ow) is observed 
to be approximately 0.0. 



4-3.1 Comparison With Other Studies 

The stepwise LF estimates for each field are mutually con- 
sistent within their uncertainties. Further, estimates from 
HUDF appear consistent with those of other authors (Oesch 
et al. 2010 and McLure et al. 2010, with only Oesch et al. 
2010 shown in Figure [T3] for clarity) while the ERS con- 
straints are consistent with the brighter sample of Ouchi et 
al. (2009) obtained using Subaru/Suprime-Cam. 

In order to explore evolution in the UV luminosity func- 
tion, we compare our results with those at lower-redshift 
using the Schechter (1976) function fits to the (z) = 3.05 
(Reddy & Steidel 2009) and (2) = 6.0 (Bouwens et al. 2006) 
luminosity functions (included in Figure I13|l . The position 
of these LFs with respect to our estimates indicates a clear 
decline in the UV LF from 2 = 3^7, together with signif- 
icant evidence of evolution from 2 = 6^7 (given that our 
stepwise LF estimates all lie below that of the Bouwens et 
al. 2006 (2) = 6.0 LF). 



4-3.2 Schechter Function Parameterisation 

To further explore the implications of the observed 2 ~ 7 lu- 
minosity function we fit the stepwise LF using the Schechter 
(1976) functional form of the luminosity function, 



(2) 



Unfortunately our stepwise LF estimates fail to extend suffi- 
ciently at lower luminosities to constrain the faint end slope 
of Schechter function. With this in mind we determine the 
remaining parameters, Af^goo and 0*, assuming a selection 
of different faint-end slope values consistent with lower red- 
shift determinations a (s {1.5, 1.7, 1.9}). The best fit values 
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M*goo r a P1600 [1025 ergs-1 Mpc'S Hz"!] (p. [Mq yr"! Mpc-^]) 

A/1600 < -18.5 (SFR > 1.5M0yr-i) < -13.0 (> 0.01 Mgyr"!) < -8.0 (> 0.0001 Mq yr"!) 

-19.8 0.00126 -1.5 3.6 ( 0.005 ) 7.8 ( 0.010 ) 8.2 ( 0.011 ) 

-19.9 0.00106 -1.7 3.6 ( 0.005 ) 10.6 ( 0.014 ) 12.7 ( 0.016 ) 

-20.1 0.00072 -1.9 3.6 ( 0.005 ) 14.9 ( 0.019 ) 25.1 ( 0.033 ) 

Table 6. The best fit values of A^jgQQ and 0» assuming fixed a G {1.5, 1.7, 1.9} together with the UV luminosity and star formation rate 
densities determined by integrating down to various depths. 
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Figure 13. The z 7 rest-frame UV (1600A) LF derived from 
each of the WFC3 fields (blue/purple points) together with con- 
temporary and lower-redshift comparisons. The solid blue line 
denotes the best fit Schochter function assuming a = —1.7. Solid 
grey lines denote the luminosity function at [z) = 2.3 (Reddy & 
Steidel 2009) and (2:) = 6.0 (Bouwens et al.). The grey circles are 
the estimates of Oesch et al. (2010) from their HUDF sample and 
from the wide area sample Ouchi et al. (2009). 
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Figure 14. 68% and 95% confidence contours for Af^ggg and 0» 
for each of the three assumed values of a. The crosses mark the 
best- fit parameters found a,t z ^ 3.05 by Reddy & Steidel (2009). 



of MigQQ and (j)* assuming each of these values of a are 
listed in Table [6] and the 68% and 95% confidence intervals 
are shown in figure [141 
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Figure 15. The UV luminosity density (left axis) and star forma- 
tion rate density (right axis) as a function of the Migoo magnitude 
down to which the luminosity function is integrated for the best 
fit luminosity functions assuming a £ {1.5, 1.7, 1.9}. The shaded 
grey box denotes the observed region with the remainder inferred 
from extrapolation of the LF. 



4-3.3 UV Luminosity and Star Formation Rate Density 

An important quantity derived from the luminosity func- 
tion is the ultraviolet luminosity density pieoo = / L x 
4'i60o{L) dL, from which the ionising photon luminosity and 
the star formation rate density can be inferred under a set 
of assumptions. Figure [15] shows the UV luminosity density 
as a function of the magnitude down to which the luminos- 
ity function is integrated for each of the best-fit luminosity 
functions with fixed a G {1.5, 1.7, 1.9}. Table [6] includes the 
luminosity density for each of these LFs for Mieoo < —18.5 
(i.e. the observed luminosity density), M1600 < —13.0 and 
A/1600 < -8.0. 

Ultraviolet emission from star forming galaxies is dom- 
inated by massive (m > 5Mq) short-lived stars. The short 
main-sequence lifetimes of these stars implies their presence, 
and thus UV emission, is a good tracer of the amount of on- 
going star formation. The star formation rate of a galaxy is 
often related to the UV luminosity through a linear equa- 
tion, i.e. 



SFR^Bvv X 10"^**!,, [ergs"' Hz"'], 



(3) 



where the parameter Buv is a calibration factor dependent 
on the previous star formation history, IMF and metallicity. 

Popular choices for this calibration parameter include 
Buv = 1.4 (Kennicutt 1998) and Buv = 1.25 (Madau, 
Pozetti & Dickinson 1998). Using the Starburst99 (Lei- 
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therer et al. 1999) population synthesis model assuming 
solar metallicity, a Salpeter IMF and 100 Myr of previous 
star formation we obtain a value of Buv = 1.31, similar 
to both the Kennicutt (1998) and Madau, Pozetti & Dick- 
inson (1998) calibrations. Applying this calibration factor 
we convert the UV luminosity density assuming the three 
luminosity functions obtained in the previous section to de- 
termine the star formation rate density. The relationship 
between star formation rate density and the magnitude to 
which the LF is integrated is shown in Figure [15] (on the 
right axis). Table [6] also includes the star formation rate 
densities (where SFR > G {1.5, 0.01, 0.0001} M© yr~^) for 
each of three best-fit luminosity function. 

Using estimates of the UV luminosity function at lower 
redshifts (z G {3.8,5.0,5.9}, from Bouwens et al. 2007) 
we can determine the dust-uncorrected high-redshift cosmic 
star formation rate density history. This is shown in Figure 
1161 again for 3 different i\/i6oo integration limits. Irrespec- 
tive of the integration limit there is a clear decline in the 
star formation rate density z — 4 ^ 7, assuming a < —1.7. 
For steeper faint-end slopes the trend becomes somewhat 
less pronounced. Estimates of the dust attenuation from the 
distribution of UV spectral slopes (e.g. Bouwens et al. 2009, 
Bouwens et al. 2010b, Wilkins et al. in prep) suggest de- 
creasing levels of attenuation at higher-redshift. Such a trend 
would further accentuate the declining trend in the star for- 
mation rate density over z = 4 — ^ 7. 

4-3.4 Reionisation 

As a result of the decreasing ambient temperature free elec- 
trons and protons recombined to form neutral hydrogen at 
z « 1100. Subsequently, the neutral hydrogen was reionised 
by a source (s) producing a sufficient number of ionising pho- 
tons (i.e. with A < 912A for hydrogen). The identity of 
these sources, and the exact epoch in which transition took 
place, remain some of the key outstanding questions in as- 
trophysics. One potential source of the ionising photons is 
the nascent z > 7 star forming galaxy population, where the 
ionising photon flux comes from a very-massive short-lived 
stellar population. 

The number of ionising photons can be expressed as a 
linear equation where the parameter Dion is the calibration 
relating the UV I6OOA luminosity to the ionising photon 
escape rate, 



A^ionpro — ^ioi 



X 10" L16OO. 



(4) 



Much like the UV luminosity to SFR calibration (Buv) Dion 
is senstive to a number of factors including the previous 
star formation duration, the metallicity and the initial mass 
function. Assuming the same scenario as used in deriving the 
star formation rate density (i.e. 100 Myr previous star for- 
mation, Z = 0.02, and a Salpeter IMF) the Starburst99 
population synthesis model implies a value of Dion = 1.2. 
Decreasing the previous star formation duration to 10 Myr 
and reducing the metallicity by a factor of 50 yields results 
an increase in the calibration to Dion = 3.1. Further, adopt- 
ing a radically different IMF, such as truncating the forma- 
tion of stars below 50 Mq can produce a further increase by 
a factor of x2. 

However, not all of the ionising photons are able to es- 
cape the immediate confines of the star forming region. Some 
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Figure 16. The evolution of the very high-rcdshift {z > 
3) star formation rate density assuming various values 
of the limiting magnitude/star formation rate (SFR > S 
{1.5, 0.01, 0.0001} Mgyr-i). z g {3.8,5.0,5.9} star formation 
rate densities are inferred from the Bouwens et al. (2007) UV 
luminosity functions. The z 7 points are those inferred 
from the best-fit Schechter functions listed in Table |6] For 
SFR >1.5 Mq yr~^ the star formation rate density is the same 
irrespective of the choice of faint-end slope. 



fraction (1 — /esc) are absorbed by the local ISM, leaving a 
reduced number that escape and are capable of ionising the 
surround neutral hydrogen. 



A^i 



ion CSC 



— fcsc ^ ^io 



(5) 



To consider whether the Universe is ionisised we employ 
a modified version (to compensate for our choice of cosmol- 
ogy) of equation 26 of Madau, Hardt & Rees (1999). This 
describes the number of ionising photons required to main- 
tain the ionisation of the Universe as a function of redshift 
and the hydrogen clumping factor C, 



A^io 



2.51 X 10**^C(1 



(6) 



Given that we infer Aionpro from the UV luminosity 
density around I6OOA and have limited constraints on both 
C and /esc it is convenient to ascertain the range of val- 
ues of C / fcsc for which the Universe can remain ionised 
assuming the luminosity functions for 3 different faint end 
slopes obtained in Section 14.31 Figure [17] shows the limit 
down to which the various luminosity functions need to be 
integrated to sustain the ionisation as a function of the 
parameter C//csc. Assuming a lower integration limit of 
Afieoo = —18.5, i.e. including only the contribution from 
galaxies brighter than our observed limit, C / fcsc is found 
be < 5.0. Assuming a faint-end slope of a = —1.7 and ex- 



14 Stephen M. Wilkins, et al. 




J I I I I I I I I I I I I I 

10 20 30 



Figure 17. The value of A/isoo it is necessary to integrate the 
best fit UV luminosity functions assuming fixed a S {1.5, 1.7, 1.9} 
to maintain the ionisation of the Universe as a function of C/fesc, 
i.e. the ratio of the clumping factor to the escape fraction of ion- 
ising photons. 

trapolating to Mieoo < —8.0 limits the value of C//csc to 
< 17.0. 

While the values of C and /esc remain uncertain a num- 
ber of studies have attempted to constrain them through 
both observations and simulations. Early simulations sug- 
gested a value of C « 30 (Gnedin & Ostriker 1997), while 
more recent attempts have yielded considerably smaller val- 
ues (e.g. Pawlik et al. 2009: C ~ 5). The value of /osc is sim- 
ilarly poorly constrained; recent observations (e.g. Siana et 
al. 2010) at z « 1-2 suggest a values of 0.05-0.20, while 
simulations (at higher-redshift) suggest potentially larger es- 
cape fractions. Assuming /esc = 0.2, and C = 5 suggests 
C//csc = 25. This value is in excess of that capable for 
the directly observed population of galaxies (i.e. those with 
Afieoo < —18.5) suggesting that we have not observed the 
full population responsible for reionisation. However, it is 
only mildly in excess of that possible if the luminosity func- 
tion is extrapolated to fainter magnitudes assuming a faint- 
end slope of Of = —1.7. Slopes steeper than a — —1.86 would, 
assuming there was a contribution from galaxies down to 
Afieoo ~ —8.0, provide sufficent ionising photons to main- 
tain the ionisation of the Universe ai z ~ 7. Of course, 
reducing the metallicity or changing the IMF, would, for a 
given UV luminosity density increase the number of ionising 
photons. Thus, at 2; « 7, assuming /esc > 0.2, and C = 5 it 
appears that star forming galaxies could produce a sufficient 
number of ionising photons assuming a fairly steep faint-end 
slope and/or lower metallicity or different IMF. 



5 CONCLUSIONS 

In this work we have presented a catalogue of candidate z- 
drop z ^ 7 galaxies from 4 separate fields (HUDF, P34, P12, 



ERS) covering a total area of ~ 50.0 arcmin^. The deepest 
field, the HUDF, extends to Ji25w ~ 28.36 {7a) and shal- 
lowest, the ERS, Ji25w = 27.03 (7a). In total we present 44 
candidates (HUDFTl, P34:15, P12:7, ERS:11) covering an 
apparent magnitude range of Ji25w = 25.9 — 28.36. Simula- 
tions reveal these galaxies likely cover a redshift distribution 
z = 6.4 — 7.7. Of these 44 candidates, 15 are recorded in a 
range of previous studies. 

While our selection criteria guards against any intrinsic 
interloper population it is still susceptible to contamination 
through photometric scatter and transient phenomena (su- 
pernovae and significant apparent motion objects). Simula- 
tions testing the effect of photometric scatter suggest a con- 
tamination rate of ~ 5% (~ 2.5 objects). Further, we note 
the presence of a possible supernovae in the HUDF field and 
record the presence of an object, which due to significant 
proper motion, mimics the characteristics of a high-redshift 
galaxy. 

Using this new data set, we determine the stepwise rest- 
frame ultraviolet luminosity function for each field (—21.0 < 
Miaoo < —18.5 over the 4 fields). These estimations of the 
luminosity function are consistent with previous estimates 
at 2 ~ 7. Further, these appear to confirm that the observed 
decline in the bright-end of the UV luminosity function (and 
star formation rate density) with redshift seen between z = 3 
and z = 6 continues to z « 7. 

The current data does not extend to sufficiently low 
luminosities to probe the faint-end slope of the luminosity 
function (assuming the Schechter parameterisation) . How- 
ever, by assuming various faint-end slopes —1.7 > a > —1.9, 
motivated by observations at lower redshift, the remaining 
luminosity function parameters can be determined. From 
these the ultraviolet luminosity density, and thus both the 
star formation rate and ionising photon luminosity densities 
can be estimated. 

The ionising photon luminosity density inferred from 
the observed luminosity function extrapolated to Afieoo = 
—8.0 assuming either a steep steep faint-end a < —1.86 or 
a shallower slope with a lower metallicity and/or different 
IMF is sufficient to maintain the ionisation of the Universe 
for recent estimates of the hydrogen clumping factor and 
ionising photon escape fraction. 
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